bcd/^ot  "^UMENTATION  PAGE 

Mk  a*  s*  \  nt.w*tN  ’o  '  ’ou/  '*w« 

■  lw  f\  "I  #  U  J  ^  tt»^9  »*c  »«■•  •**•*9  :*#  .  ;neci  on  of  .Mo'*««t» 

^<«aM  #  W<|«9  tf»»*  Dura#**  :  2V4»h.rqton  M«*OQu4rt*’ 


fo/’m  Approved 
OMB  No  0704  01S8 


C1*  %  *0  9««M9#  !  'o«f  +\OOr\4  nc'u<3iP9  t*t  itm*  fo*  f*v>*»Mr»g  '"%ir  ••  \t  *»  ;  j»?  «  -,c 

tt»^9  »««  »*.-e*Mr»q  .enter  on  of  in«o'***«tion  S#n<|  commtmt  reajramg  thuOw'd**  estimate  5/  <n«  ;»n#r  j-»&«u 
uonq  lf»*s  o^ra«n  !;  AfAsnrpqtori  ***0<|u4r|*r*  S*rv>ctv  O'ftCtOrttff  for  inforff**1>on  Oo^r«t<on>  <r.<j  o*oort%  ' i  S  «•*♦ 
*nq  1C  t*»e  0**  •*  '  4f»q  8ud<J*l  »4pef^on»  H*durtiO«P'O|«<U0fC4  01»IJ  Aw.nqtcn  1C  2CS0I 


2  REPORT  DATE 

April  22,  1994 


4.  IIILt  ANU  3UOIIIU 


Anionic  and  Neutral  States  of  Li30 


6.  AUTHOR(S) 

Maciej  Gutowski  and  Jack  Simons 


7.  PERFORMING  ORGANIZATION  NAME(S)  ANO  ADORE SS(ES) 

UNIVERSITY  OF  UTAH 
DEPARTMENT  OF  CHEMISTRY 
SALT  LAKE  CITY,  UTAH  84112 
U.S.A. 


9.  SPONSORING  /MONITORING  AGENCY  NAME(S)  ANO  AOOR 

OFFICE  OF  NAVAL  RESEARCH 
CHEMISTRY  PROGRAM 
800  NORTH  QUINCY  ST. 

ARLINGTON,  VIRGINIA  22217-5000 


11.  SUPPLEMENTARY  NOTES 


3.  REPORT  TYPE  ANO  OATES  COVERED 

TECHNICAL  REPORT 


S.  FUNDING  NUMBERS 

G  N00014-89-J-1497 
R&T  4131050 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


Technical  Report  No.  45 


10.  SPONSORING /MONITORING 
AGENCY  REPORT  NUMBER 


Submitted  for  publication  in  the  Journal  of  Chemical  Physics 


12a.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

THIS  DOCUMENT  HAS  BEEN  APPROVED  FOR  PUBLIC  RELEASE  AND 
SALE;  ITS  DISTRIBUTION  IS  UNLIMITED 


12b.  DISTRIBUTION  CODE 


UNLIMITED 


13.  ABSTRACT  (Maximum  200  wordi)  ...  .  . _  .  .  ,  .  .  .  .  .  , 

The  ground  state  of  the  cation,  three  lowest  electronic  states  of  the  neutral,  and  two  anionic  states  of 
Li3<D  were  studied  using  different  ab  initio  techniques.  Stationary  points  on  the  potential  energy  surfaces 
were  determined  both  at  complete  active  space  (CAS)  self-consistent  field  (SCF)  and  at  second-oider  Moller- 
Plesset  (MP2)  levels  of  theory.  Excited  states  were  approached  using  the  single-excitation  configuration 
interaction  (CIS)  method.  Electron  detachment  energies  for  the  anionic  and  neutral  states  were  calculated  at 
the  quadratic  configuration  interaction  (QCD  level  with  single,  double,  and  approximate  triple  excitations 
(SD(T))  included.  The  calculations  indicate  that  LyO*  possesses  two  bound  electronic  states.  The  ground 
1  Ai'  state  has  an  equilibrium  D3h  structure  and  a  vertical  electron  detachment  energy  (VDE)  of  0.66  eV.  The 
%'  bound  state  pseudorotates  through  3  A]  and  3B2  stationary  points.  The  barrier  for  pseudorotation  was 

found  to  be  less  than  0.002  eV  at  the  QCISDCD  level.  Two  VDE  peaks  for  the  3E'  anion  were  predicted  to  be 
at  0.45  and  1.15  eV,  for  transitions  to  the  ground  and  the  first  excited  state  of  the  neutral,  respectively.  The 
ground  state  of  the  cation  and  the  first  three  electronic  states  of  the  neutral  LijO  were  also  considered  and 

the  vertical  ionization  potential  for  the  ground  neutral  state  was  found  to  be  3.60  eV.  Li30  and  Li30"  are 
thermodynamically  stable  with  respect  to  the  unimolecular  decompositions  Li30^°^  -»  Li20  +  Li^°\  Hence 
the  species  should  be  amenable  to  experimental  studies. 


14.  SUBJECT  TERMS 

% 

17.  SECURITY  CLASSIFICATION 

OF  REPORT 

UNCLASSIFIED 

18.  SECURITY  CLASSIFICATION 
OF  THIS  PAGE 

UNCLASSIFIED 

19.  SECURITY  CLASSIFICATION 

OF  ABSTRACT 

UNCLASSIFIED 

IS.  NUMBER  OF  PAGES 


16.  PRICE  CODE 


Standard  Form  298  (Rev  2  89) 

Presented  by  Md  9 

298  102 


OFFICE  OF  NAVAL  RESEARCH 


Contract  N00014-89-J-1497 

R&T  Code  4131050 
Scientific  Officer:  Dr.  R.  Nowak 

Technical  Report  No.  45 


"Anionic  and  Neutral  States  of  U3O" 
by 


Maciej  Gutowski  and  Jack  Simons 


Prepared  for  Publication  in 


The  Journal  of  Chemical Physics 


The  University  of  Utah 
Department  of  Chemistry 
Salt  Lake  City,  Utah  841 1 2-1194 


April  22,  1994 


Accesion  For 

NTiS  CRA&I 
OTIC  TAB 
Unannounced 
Justification 


By . 

Distribution  / 


Availability  Codes 

Dist 

Bd 

Avail 

Spe 

and  /  or 
‘cial 

-  . 

1 


Reproduction  in  whole  or  in  part  is  permitted  for  any 
purpose  of  the  United  States  Government 

This  document  has  been  approved  for  public  release  and 
sale;  its  distribution  is'unlimited. 


I 

\ 


rrw< 

S-Lifiw  I  fTf  0 


Anionic  and  Neutral  States  of  Li30 

Maciej  G*itowski@  and  Jack  Simons 

Chemistry  Department 
University  of  Utah 
Sait  Lake  City,  UT  84112 


Abstract 


The  ground  state  of  the  cation,  three  lowest  electronic  states  of 
the  neutral,  and  two  anionic  states  of  Li30  were  studied  using 
different  ab  initio  techniques.  Stationary  points  on  the  potential 
energy  surfaces  were  determined  both  at  complete  active  space 
(CAS)  self-consistent  field  (SCF)  and  at  second-order  M0ller-Plesset 
(MP2)  levels  of  theory.  Excited  states  were  approached  using  the 
single-excitation  configuration  interaction  (CIS)  method.  Electron 
detachment  energies  for  the  anionic  and  neutral  states  were 
calculated  at  the  quadratic  configuration  interaction  (QCI)  level  with 
single,  double,  and  approximate  triple  excitations  (SD(T))  included. 
The  calculations  indicate  that  LijO*  possesses  two  bound  electronic 
states.  The  ground  1 A  i'  state  has  an  equilibrium  D3h  structure  and  a 
vertical  electron  detachment  energy  (VDE)  of  0.66  eV.  The  3E’  bound 
state  pseudorotates  through  3Aj  and3B2  stationary  points.  The 
barrier  for  pseudorotation  was  found  to  be  less  than  0.002  eV  at  the 
QCISD(T)  level.  Two  VDE  peaks  for  the  3E’  anion  were  predicted  to  be 
at  0.45  and  1.15  eV,  for  transitions  to  the  ground  and  the  first 
excited  state  of  the  neutral,  respectively.  The  ground  state  of  the 
cation  and  the  first  three  electronic  states  of  the  neutral  Li30  were 
also  considered  and  the  vertical  ionization  potential  for  the  ground 
neutral  state  was  found  to  be  3.60  eV.  Li30  and  Li30*  are 
thermodynamically  stable  with  respect  to  the  unimolecular 
decompositions  Li30*(°)  -»  Li20  +  LH°).  Hence  the  species  should  be 
amenable  to  experimental  studies. 


@  Present  address:  IBM  Research  Division,  Almaden  Research  Center, 
650  Harry  Road,  San  Jose,  CA  95120-6099. 
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I.  Introduction 


Recent  theoretical  studies  indicate  that  neutral  molecular 
radicals  containing  alkali  metal  atoms  can  accommodate  more  than 
one  bound  anionic  state.  The  simple  alkali  metal  oxide  diatomic 
molecules  possess  electronically  bound  anionic  states  of  3n,  1n,  1X+, 
and  3I+  symmetry  [1].  Experimental  [2,3]  and  theoretical  [4]  studies 
on  alkali  metal  trimer  anions  (M3-)  have  concentrated  primarily  on 
the  lowest  isomer.  Recently,  we  demonstrated  that  one  triplet 

(3A2‘)  and  two  quintet  (5Ai"  and  5A2')  states  of  Li3‘  and  Na3*  are  also 
electronically  stable  for  a  wide  range  of  molecular  geometries  [5]. 
Also  LiFLi-  possesses  two  electronically  bound  states  of  and  3IU+ 
symmetry  [6]. 

The  neutral  Li30  has  already  been  studied  both  experimentally 
[7]  and  theoretically  [8-11].  Wu  at  al.  identified  Li30  in  the  course  of 
their  study  on  the  nature  of  gaseous  species  over  solid  lithium  oxides 
[7].  The  values  of  the  ionization  potential,  atomization  energy,  and 
dissociation  energy  to  produce  Li  +  Li20  were  found  to  be  4.54±0.2 
eV,  228. 7±2  kcal/mol,  and  50.7±10  kcal/mol,  respectively  [7,12], 

Li30  is  an  example  of  a  "hypermetalated"  molecule  with  a 
stoichiometry  which  violates  the  octet  rule  [8].  In  addition,  it  is  a 
promising  candidate  for  being  a  "superalkali",  (i.e.,  a  molecular 
system  whose  first  ionization  potential  is  smaller  than  that  of  the  Cs 
atom  [9-11]).  Both  of  these  unusual  chemical  features  are  related  to 
the  nature  of  the  highest  occupied  molecular  orbital  (HOMO).  It 
displays  bonding  interactions  between  pairs  of  Li  "ligands"  which 
help  to  offset  the  octet-rule-violating  structure  and  antibonding  Li-0 
interactions. 

Early  predictions  of  a  Cj\  equilibrium  structure  for  Li30  [8] 
have  not  been  confirmed  in  more  advanced  calculations  [13,10] 
which  indicated  a  D3h  structure,  similar  to  that  of  Li30+.  Theoretical 
predictions  of  the  adiabatic  and  vertical  ionization  potential  produced 
3.55  and  3.45  eV,  respectively  [10,11],  far  outside  the  range  of 
experimental  values  [7],  although  the  theoretical  dissociation  energy 
of  47.1  kcal/mol  [10]  agrees  well  with  the  experimental  prediction. 


To  the  best  of  our  knowledge,  neither  the  excited  electronic 
states  of  the  neutral  nor  the  anionic  states  of  Li30  have  yet  been 
experimentally  studied. 

II.  Computational  Aspects 

For  the  lithium  atom,  we  used  the  Dunning  (9s5p/3s2p)  one- 
electron  basis  set  [14]  supplemented  with  diffuse  s  and  p  functions 
with  the  same  exponent  0.0074  [15]  and  one  d  function  with  the 
exponent  0.2  [16].  This  basis  set  is  detailed  in  Ref.  [5].  For  the  oxygen 
atom  we  employed  Dunning's  aug-cc-pVDZ  basis  set  which  was 
designed  to  describe  anionic  species  [17].  Cartesian  d  functions  were 
used  throughout  the  calculations  and  the  full  basis  set  for  the 
molecule  consists  of  82  contracted  gaussian  functions. 

Potential  energy  surfaces  were  explored  within  a  complete 
active  space  (CAS)  self-consistent  field  (SCF)  formalism  as  well  as  at 
the  second-order  Moller-Plesset  (MP2)  theory  level. 

In  CAS  SCF  calculations  we  imposed  a  constraint  that  molecular 
orbitals  which  result  from  the  core  Is  atomic  orbitals  were  doubly 
occupied  in  every  configuration  state  function  (CSF).  The  neglected 
core-core  and  core-valence  correlation  effects  are  negligible  for  the 
lithium  and  oxygen  atoms  due  to  the  low  polarizability  of  the  Is 
cores.  The  remaining  eight  (cation),  nine  (neutral)  or  ten  (anion) 
electrons  were  distributed  in  all  possible  ways  among  four  ai,  two  bi, 
and  two  b2  molecular  orbitals  (the  C2v  symmetry  labeling  is  used). 
This  choice  of  the  active  space  led  to  1764,  2352,  1176,  and  1512 
CSF's  for  the  cation,  neutral,  singlet  and  triplet  anion,  respectively. 
The  CAS  SCF  calculations  were  performed  with  the  Utah  MESS-KIT 
modular  electronic  structure  codes  [18]  which  generate  analytical 
second  geometrical  derivatives.  Stationary  points  on  the  potential 
energy  surfaces  were  determined  using  our  automated  surface 
walking  algorithm  [19]. 

In  the  case  of  the  neutral  22Ai  state,  which  correlates  with  the 
2E'  state,  the  CAS  SCF  optimization  was  hindered  by  the  problem  of 
"root-flipping"  [20].  Since  the  MP2  approach  is  inapplicable  for 
excited  electronic  states,  we  invoked  single-excitation  configuration 
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interaction  (CIS)  approach  [21]  to  determine  geometry  and  relative 
energy  of  the  22Ai  transition  state.  The  CIS  results  for  the  doublet 
neutral  state  should  be  considered  cautiously.  In  the  CIS  approach, 
creation  of  spin  eigenstates  relies  on  having  an  RHF  ground  state  and 
non-interaction  of  singlets  and  triplets  does  not  carry  over  to 
doublets  and  quartets  with  a  UHF  reference.  The  Gaussian  92  code 
[22]  attempts  to  handle  these  doublet  cases  but  the  theory  is  not 
clean  anymore  [23].  In  view  of  the  above  doubts,  we  "calibrated"  the 
CIS  approach  on  the  2B2  electronic  state  which  also  correlates  with  2E' 
but  for  which  the  CAS  SCF  and  MP2  approaches  are  straightforwardly 
applicable. 

In  the  MP2  geometry  optimizations  we  allowed  for  the 
correlation  of  the  core  orbitals.  In  general,  the  structures 
corresponding  to  stationary  points  are  quite  similar  at  the  CAS  SCF 
and  MP2  levels  what  suggests  that  the  core-core  and  core-valence 
correlation  effects  are  not  important  for  geometrical  predictions. 

The  restricted  CAS  SCF  approach  is  capable  to  predict  accurate 
geometries  but  it  is  inappropriate  to  accurately  compare  energies  of 
species  with  a  different  number  of  electrons.  Hence  we  employed  the 
quadratic  configuration  interaction  (QCI)  approach  with  single, 
double  and  approximate  triple  excitations  (SD(T))  [24]  to  determine 
relative  energies  and  electron  detachment  energies.  The  QCISD(T) 
approach  is  size-extensive  and  takes  into  account  dynamical 
correlation  affects.  In  the  QCI  calculations,  we  kept  the  eight  core 
electrons  uncorrelated.  We  checked  in  our  earlier  study  on  LiFLi*  [6] 
that  such  a  restriction  changes  the  vertical  detachment  energy  for 
the  ground  state  anion  by  less  than  0.003  eV.  The  QCI  results  were 
obtained  with  the  Gaussian  92  suite  of  codes  [22]. 

III.  Results 

The  stationary  points  on  the  potential  energy  surfaces  of  the 
cation,  neutral,  and  anion  determined  at  -the  CAS  SCF  and  MP2  levels 
are  characterized  at  Table  1  and  the  geometrical  parameters  used  are 
defined  in  Fig.  1.  The  vibrational  frequencies  were  calculated  using 
analytical  second  derivatives.  For  the  excited  2B2and  22Ai  states  of 


the  neutral,  stationary  points  are  also  reported  at  the  CIS  level.  The 
spatial  extents  of  the  electronic  charge  distributions  are 
characterized  by  the  SCF  values  of  <R2>.  The  reported  relative 
energies  were  obtained  at  the  QCISD(T)  level,  as  were  our  vertical 
(VDE)  and  adiabatic  electron  affinities  (EAa)  presented  in  Tables  2 
and  3  for  the  neutral  and  the  anion,  respectively. 

A.  The  Cation  and  Neutral 

The  D3h  electronic  configuration  for  the  closed-shell  cation 
involves  an  orbital  occupation  3ai’22e'2la2"2  (the  core  Is  orbitals  are 
not  included  in  this  labelling)  which  will  be  denoted  (+).  Our  D3  h 
geometry  and  frequencies  for  the  closed-shell  cation  are  quite 
similar  at  the  CAS  SCF  and  MP2  levels.  They  are  also  in  a  good 
agreement  with  the  results  of  Rehm  et  al.  [10].  The  piramidization 
mode  (a2"),  which  lowers  symmetry  to  C3V,  has  the  lowest  frequency. 
Interestingly,  the  cation  H3  0+,  studied  earlier  [25]  has  a  C3  v 
equilibrium  structure  with  the  D3h  barrier  of  only  1.07  kcal/mol  [26]. 

The  MP2  equilibrium  structure  of  the  ground  state  neutral  with 
the  electronic  configuration  (+)4ai’  is  again  D3h  and  similar  to  that  of 
the,  cation.  Because  the  2Ai'  wavefunction  is  found  to  have  a  single 
configuration  character,  the  MP2  prediction  should  be  quite  accurate. 
Unexpectedly,  the  e’  (in-plane)  vibrational  mode  has  an  imaginary 
frequency  at  the  CAS  SCF  level.  Doubting  whether  the  CAS  SCF  e’ 
imaginary  frequency  is  physically  meaningful,  we  carried  out  a  CAS 
SCF  search  for  a  new  stationary  point  in  C2v  symmetry  and  found  a 
structure  with  Ri  =  1 .683  A,  R2=1.731  A,  0=133°  and  an  energy  0.06 
eV  lower  than  at  the  D3h  stationary  point.  However,  the  QCISD(T) 
energy  is  lower  by  0.02  eV  at  the  D3h  than  at  the  C2v  CAS  SCF 
stationary  point.  Hence,  we  conclude  that  the  D3h  structure 
corresponds  to  a  genuine  minimum  and  the  CAS  SCF  approach  suffers 
for  symmetry  breaking  artifacts  [27].  Even  though  this  2Ai'  state  has 
a  D3h  equilibrium  structure,  the  MP2  -vibrational  frequencies  of  the 
a2"  and  e'  modes  are  much  softer  than  in  the  underlying  cation  which 
implies  that  C3V,  C2v,  and  Cs  symmetries  are  easily  accessible  by 
vibrational  movement. 


6 


The  4aj'  HOMO  of  the  neutral  is  dominated  by  Li  2s  orbitals 
which  interact  constructively  with  each  other  and  destructively  with 
small  s-type  contributions  from  the  central  O  atom,  as  observed  in 
Refs.  [8-10].  The  Mulliken-population  atomic  charges  are  +0.20  and 
-0.60  for  Li  and  O,  respectively.  Our  values  for  the  VDE  and  EAa  of 
this  ground  state  of  Li30  are  3.60  and  3.59  eV,  respectively,  in  a  good 
agreement  with  the  electron  propagator  theory  VDE  of  3.45  eV  [11]. 
However,  these  results  disagree  with  the  experimental  estimation  of 
4.54±0.2  eV  [7].  The  theoretical  VDE  for  Li30  is  lower  than  that  of 
alkali  metal  atoms,  so  our  result  support  the  claimed  "superalkali" 
nature  of  Li30  [9-11]. 

The  first  excited  state  for  the  neutral  of  2E'  symmetry  has  a 
dominant  orbital  occupancy  of  (+)3e'  and  is  subject  to  first-order 
Jahn-Teller  (FOJT)  distortion.  Geometry  optimization  for  the  resulting 
2B2  component  of  this  state  is  straightforward  and  the  resulting  CAS 
SCF  and  MP2  geometries  and  frequencies  are  quite  similar.  Due  to 
problems  with  "root-flipping"  at  the  CAS  SCF  level,  geometry 
optimization  for  the  22Ai  component  of  this  2E’  state  was  performed 
at  the  CIS  level.  The  reliability  of  the  CIS  approach  was  tested  on  the 
2B2  state  where  the  stationary  point  characteristics  were  found  to  be 
similar  for  the  CIS,  CAS  SCF,  and  MP2  approaches.  We  therefore 
believe  that  application  of  the  CIS  approach  to  the  22A  i  state  is 
justified. 

Pseudorotation  in  the  2E'  state  is  depicted  in  Fig.  2.  Estimating 
the  pseudorotation  barrier  requires  a  consistent  calculation  of  the 
22Ai  and  2B2  energies.  A  QCISD(T)  estimate  for  the  22Ai  energy  was 
obtained  using  the  following  approximation:  Eqci(2B2)  -  Ecis(2B2)  + 
Ecis(22Ai),  because  only  the  CIS  method  could  be  applied  directly  to 
the  22A  i  state.  The  barrier  for  pseudorotation  thus  found  was  less 
than  0.04  eV. 

The  CIS  oscillator  strength  for  the  2E’  «-  2Aj'  transition  is  0.48 
and  the  QCISD(T)  vertical  excitation  energy  is  0.728  eV.  Interestingly, 
the  neutral  2E'  state  thus  produced  would  be  in  the  neighborhood  of 
its  conical  intersection.  This  fact  could  be  reflected  in  the  dynamics  of 
the  2E’  state,  which  might  be  studied  using  time  dependent  two- 
photon  ionization  techniques  [28]. 
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The  next  excited  state  of  L13O  has  2 A 2"  symmetry  and  (+)2a2M 
orbital  occupancy  and  an  equilibrium  D3h  structure.  In  comparison 
with  the  cation,  the  vibrational  frequencies  are  quite  similar  and  the 
equilibrium  R  is  somewhat  shorter.  These  features  are  consistent 
with  the  "out-of-the-plane"  nature  of  the  2a2"  orbital  and  similar 
trends  were  observed  in  the  alkali  metal  trimers  with  the  unpaired 
electron  in  an  a2"  orbital  [5,29].  In  Li30,  the  2a2"  orbital  is  dominated 
by  the  Li  2p  atomic  orbitals  with  constructive  Li-Li  and  destructive 
Li-O  interactions.  The  CIS  oscillator  strength  for  the  2A2"  <-  2Aj' 
transition  is  0.32  and  the  QCISD(T)  vertical  excitation  energy  is  1.551 
eV. 

B.  The  Anion 

The  ground  electronic  state  of  Li30*  has  a  D3h  equilibrium 
structure  with  a  geometry  close  to  that  of  the  cation  and  that  of  the 
ground  state  neutral.  There  is  again  a  significant  difference  between 
the  a2"  out-of-plane  vibrational  frequency  calculated  at  the  MP2  and 
CAS  SCF  levels,  but  this  time  both  approaches  predict  a  D3  h 
minimum.  The  MP2  frequencies  are  quite  similar  for  the  anion  and 
the  2Ai'  neutral. 

The  wavefunction  for  the  anionic  1 A 1'  state  is  dominated  by  the 
(  +  )4ai'2  configuration,  but  two  equivalent  contributions  from  the 
(+)3e'2  configuration  are  also  important  and  have  Cl  coefficients  of 
0.30  each.  This  feature  is  consistent  with  the  low  separation  between 
the  2 A 1 '  and  2E’  states  of  the  neutral  (0.7  eV)  and  suggests  that  3E’ 
anionic  state  may  be  also  electronically  stable.  The  VDE  for  the  1 A 1' 
anion  state  is  0.656  eV  and  the  electron  detachment  peak  is  expected 
to  be  sharp  because  the  anion  and  neutral  equilibrium  structures  are 
very  similar. 

In  addition  to  the  ground  1 A 1'  state,  the  anion  Li30*  possesses  a 
second  electronically  bound  state  of  3E’  symmetry  with  the  dominant 
electronic  configuration  (+)4ai'3e\  Due-  to  FOJT  distortion,  stationary 
points  develop  on  3B2  and  3Ai  surfaces.  Both  CAS  SCF  and  MP2 
approaches  locate  a  transition  state  on  the  3B2  surface  with  negative 
curvature  along  the  b2  distortion  mode.  The  surface  must,  however. 
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be  extremely  flat  since  at  the  QCISD(T)  level  the  relative  order  of  the 
3B2  and  3Ai  stationary  points  is  changed,  giving  3B2  lower  than  3Ai 
by  0.001  eV,  whereas  the  difference  at  the  MP2  level  was  0.011  eV 
with  the  opposite  order.  Clearly,  the  pseudorotation,  as  depicted  on 
Fig.  2,  is  practically  free  and  the  numerical  values  of  the  b2  mode 
frequency  are  probably  of  little  reliability. 

The  electron  detachment  energies  for  the  pseudorotating  3E’ 
state  were  calculated  at  both  daughter-state  stationary  points  (see 
Table  3).  The  detachment  energies  to  the  ground  2Aj'  state  of  the 
neutral  are  predicted  to  lie  in  the  range  0.43-0.46  eV.  The 

detachment  energies  to  the  higher  pseudorotating  2E'  state  of  the 
neutral  lie  within  1.1 -1.2  eV.  We  thus  conclude  that  the  electron 

detachment  peaks  from  the  3E'  anion  state  would  bracket  the  anion’s 

ground  state  detachment  peak  at  0.66  eV.  Hence,  the  3E’  state  may  be 

amenable  to  experimental  detection,  providing  a  significant 
concentration  of  the  triplet  can  be  produced  in  the  source. 

Another  feature  which  makes  Li30  and  LisO*  suitable  for 
experimental  studies  is  their  thermodynamic  stability.  The  energy 
barrier  (corrected  for  zero-point  vibrations)  for  the  decomposition 
Li30(+°->  U2O  +  Li(+0')  is  predicted  to  be  84.40,  44.54,  and  46.03 
kcal/mol  for  the  cation,  neutral,  and  anion,  respectively.  The 
increased  stability  of  the  anion  compared  to  that  of  the  neutral 
reflects  the  fact  that  the  electron  affinity  is  larger  for  Li30  than  for 
Li.  Our  decomposition  energy  for  the  neutral  is  within  the  range  of 
experimental  data  50.7±10  kcal/mol  [7,12]. 

IV.  Conclusions 

Theoretical  calculations  ind;cate  that  LUO*  can  possess  more 
than  one  bound  electronic  state.  The  fully  symmetric  singlet  state  is 
the  ground  state,  but  a  triplet  state  is  also  electronically  bound. 
Although  only  the  anionic  states  of  Li30*  were  studied  in  this  work, 
the  preliminary  results  for  the  isoelectroi.ic  LUN*  produce  a  VDE 
from  the  *Ai  state  of  ca.  0.51  eV.  For  the  3T2  state  of  LUN*,  the  VDE 
is  ca.  0.17  eV,  whereas  detachment  to  the  2T2  state  of  the  neutral 
LUN  would  require  ca.  0.98  eV. 
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For  the  ground  1 A  j '  state  of  Li30*,  the  electron  detachment 
energy  is  0.66  eV.  The  cation,  ground  state  neutral,  and  the  anion  all 
have  the  equilibrium  D3h  geometries,  with  the  O-Li  distances  in  the 
range  1.71-1.69  A  (MP2  level). 

The  pseudorotating  anionic  3E'  state  is  electronically  stable  with 
respect  to  the  2 At'  and  2E’  states  of  the  neutral  by  ca.  0.4  and  1.1  eV, 
respectively  The  geometrical  features  of  the  daughter-state  3Aj  and 
3B 2  stationary  points  are  consistent  with  the  bonding/antibonding 
interactions  among  the  atomic  orbitals  contributing  to  the  anion's 
singly  occupied  molecular  orbitals.  The  QCISD(T)  energies  of  the  two 
stationary  points  are  the  same  to  within  0.002  eV.  Hence,  the 
pseudorotation  is  practically  free. 

For  the  neutral  Li30,  we  studied  the  2E’  and  2A2"  excited  states 
in  addition  to  the  -Af  ground  state.  The  2E’  state  pseudorotates 
through  a  2B2  minimum  and  a  22A  i  transition  state  with  a 
pseudorotation  barrier  of  ca  0.04  eV.  The  oscillator  strength  for  the 
2E'  <r-  2Aj’  transition  is  0.48  and  the  corresponding  vertical  excitation 
energy  is  0.728  eV. 

The  2 A2"  neutral,  with  a  D3h  equilibrium  structure,  is  separated 
from  the  ground  state  neutral  by  1.551  eV,  and  is  easily  accessible 
since  the  oscillator  strength  is  0.32. 

Our  vertical  and  adiabatic  ionization  potentials  for  Li30  of  3.60 
and  3.59  eV,  respectively,  disagree  with  the  experimental  value  of 
4.54±0.2  eV  [7],  yet  agree  with  other  theoretical  predictions  [10,11]. 

All  of  the  species  discussed  in  this  study  are 
thermodynamically  stable  with  respect  to  unimolecular 
decomposition.  Hence,  the  neutral  and  the  anion  are  more  amenable 
to  experimental  studies  than  their  hydrogen-substituted  analogs  H3O 
and  H3O',  which  are  thermodynamically  unstable  [25], 
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Captions  for  Figures 


Figure  1.  Geometrical  parameters  for  the  C2v  structure  of  the  L13O 
neutral  and  anionic  species. 

Figure  2.  Graph  representing  the  pseudorotation  for  anionic  or 
neutral  Li30. 


Ic  I.  Stationary  points  (distances  R  in  A.  i')  in  degrees)  and  harmonic  frequencies  (cm'1)  for  different  electronic  states  of 
cation,  neutral,  and  anion  of  U3O.  The  spatial  extent  of  the  SCF  electronic  charge  distribution  <R^>  is  given  in  a.u.,  and 
relative  energies  (QCISD(T)  level)  are  in  cV.  The  dominant  electronic  configurations  (DEC)  are  given  with  respect  to  the 
cd-shell  cationic  core  3ai’^2e'2  Ia2"^.  which  is  denoted  (+). 
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“'IR  intensities  (km/mol)  in  parenthesis 

b^thc  energy  of  22A|  is  estimated  as:  Eqc!(^B2)  +  EciS(2^Ai)  -  EciS(^B2) 
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Table  2.  Electron  vertical  detachment  energies  (VDE)  and  adiabatic 
electron  affinities  (EAa)  (in  eV)  calculated  at  the  QCISD(T)  level  for 
the  transitions  U3O+  +  e  <-  U3O. 


» 


Transition 

VDE 

EAa  . 

lAi'  +  e  «-  2Aj' 

3.603 

3.593 

jAi  +  e  <-  2B2 

2.933 

2.890 

*Ai  +  e  <-  22Aja) 

2.954 

2.909 

iAi’  +  e<-2A2" 

2.052 

2.042 

a)the  energy  of  22Ai  is  estimated  as  in  Table  1 


Table  3.  Electron  vertical  detachment  energies  (VDE)  and  adiabatic 
electron  affinities  (EAa)  (in  eV)  calculated  at  the  QCISD(T)  level  for 
the  transitions  U3O  +  e  <-  U3CK 


Transition 

VDE 

EAa 

2Ai'  +  e  «-  lAi' 

0.656 

0.656 

12Ai  +  e  *—  3Ai 

0.450 

0.430 

22Ai  +  e  < —  ^ a  1 

1.346a) 

1.074b) 

12Ai  +  e  <—  3B2 

0.459 

0.431 

2B2  +  e  <-  3B2 

1.151 

1.134 

a)the  VDE  is  estimated  as:  Eqci(3Ai)  -  Eqci(2B2)  +  Ecis(2B2  )  -  Ecis(22Ai). 


b)the  EAa  is  estimated  as:  VDE  +  Ecis(22Ai  at  the  3Ai  geometry)  -  Ecis(22Ai 
at  its  CIS  stationary  point). 


